Surprising properties of water on its binodal as the reflection of 
the specificity of intermolecular interactions 

V.L. KulinskiiQ 

Department of Molecular Physics, Kiev National 
University, Glushko 6, 65026 Kiev, Ukrain^ 

N.P. MalomuzhI 

Department for Theoretical Physics, Odessa National University, 
Dvoryanskaya 2, 65026 Odessa, Ukraine 

Abstract 

In the paper the behavior of density (or specific volume), the heat of evaporation and entropy 
per molecule for normal and heavy water on their coexistence curves is discussed. The special 
attention is paid on the physical nature of the similarity in the behavior of density and the heat of 
evaporation for water and argon as well as the nearest water homolog H2S. It is shown that the 
appearance of this similarity is a consequence of the rotational motion of water molecules, which 
averages the inter-particle potential in water and leads it to the argon-like form. To describe the 
fine distinctions in the behavior of the binodals for water, H2S and argon the dependence of the 
proper molecular volume on pressure is taken into account. In accordance with this often used 
the van der Waals and Carnahan-Starling equations of states are modified. The very surprising 
behavior of the entropy diameter for water is analyzed. It is shown that the nontrivial details of 
the temperature dependence for the entropy diameter is directly connected with the peculiarities 
of the rotational motion of molecules in water. The effect of strong dimerization of water molecules 
in the fluctuation region near the critical point is studied in details. 
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I. INTRODUCTION 



The surprising properties of water, as it is well known |T], are caused by H-bonds. The 
character of their spatial organization and time dynamics depends essentially on temperature 
and pressure. As a result, the peculiarities of the translational and rotational motion of water 
molecules are very different near the triple and critical points. There are several key facts 
which allow to make the important conclusions about the behavior of the H-bond network 
in water and the character of the thermal motion of water molecules, in particular, near its 
critical point. One of them is presented in Fig. where the comparison of the specific volume 
per molecule v^^^ for water and argon (i = w; Ar) is given. The dimensionless temperature 
t = T/Tc{i) , where Tc(z) is the critical temperature of liquids, is used. As seen, the behavior 
of the specific volume of water is very surprising. Practically in the whole region it is argon- 
like. Only in the narrow vicinity of the critical point (0.95 < t < 1) the deviation of v^'^\t) 
from the argon-like dependence is essential. Outside of this region (t^ < t < 0.95), where 
= 0.42 is the melting point, the deviation from the argon-like dependence does not exceed 
3 ^ 4% . The minimum of the specific volume of water at T ^ 2HK corresponds to the 
weak increase of the ratio Rv{t) = v^^\t) / v^'^^\t) at approaching the melting point. In 
accordance with Fig. || the temperature dependence of Rv(t) can be represented in the form: 



where rnit) is the contribution caused by H-bonds in water: rnit) <C k. 

To describe the temperature dependence of rnit) seems to be natural to apply to the 
Hubert's principle [? ? ], which was formulated for the first time in the algebraic invariant 
theory and having the numerous applications in the statistical hydrodynamic ^ . According 
to this principle, an arbitrary complicated function can be expanded in the series with 
respect to independent primitive functions, which have the same properties of symmetry: 



where are the coefficients, which are functions of temperature t and the dimensionless 
pressure p = P/Pc, Pc is the critical pressure. The independent structural characteristics of 
the H-bond network, so-called structural functions [? ], play a role of primitive functions 
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FIG. 1: Temperature dependence of the ratio Ry = v^'^\t)/v^^'^\t) on the coexistence curves of 
water and argon according to [? ]. 

The most important characteristics of the H-bond network are the average number nn 
of H-bonds per molecule and the parameter of tetrahedricity x(T) (see P, 0). Therefore, 
the first two structural functions can be chosen as: Si = nH{t,p) and 5*2 = XH{t,p),- The 
structural functions of a higher order are assumed to be responsible for the finer details of 
the H-bond network and here they will be ignored. In [? ] it was shown, that with good 
accuracy (|l]) can approximated as 

rH{t)=rP -nnit), n^^(t) = 4(1 - At + ...). (2) 

Fitting the experimental data on Rv{t) with the help of (0) one can get: k = 0.63, r^^ = 
0.015; A =0.85. The value and the temperature dependence of nnit), given by Eq. (^, are 
in good agreement with the results obtained in [P, P| pi|. 

Below we will show, that the temperature dependence of the heat of evaporation in water 
is also similar to that for argon. The estimate for nH{t), following from the analysis of the 
heat of evaporation, practically coincides with Eq. (H). 

To explain the similarity in the behavior of the specific volume and the heat of evapora- 
tion in water and argon we should conclude that the mentioned quantities and some other 
thermodynamic characteristics of water are determined by the averaged intermolecular po- 
tential, which takes the argon-like behavior. 
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From Eq. (0) it follows that near the critical point the average number of H-bonds per 
molecule takes value n//(l) ^ 0.7. It means, that in the fluctuation region, 0.85 < t < 1, 
water can be considered as the ensemble of dimers. This conjecture is supported by the 
results obtained in [? ]. Besides, it becomes to be self-evident if we take into account that in 
this region the volume of cavity {vcav = 2fm) occupied by two water molecules are greater 
than that volume, which is created by a rotating dimer {vcav > Vd) ■ Qualitatively, almost full 
dimerization of water molecules allows us to explain naturally the sharp enough increment 
of the ratio near the critical point (see Fig. [l|). 

The present paper is devoted to the careful analysis of the density, entropy and the heat 
of evaporation on the coexistence curve of water in the wide temperature range from the 
melting to the critical points. The special attention is paid to the dimerization effects of 
water molecules in the fluctuation region. 

To describe the fine details in the behavior of the density, diameter of entropy and the 
heat of evaporation on the binodal we take into account the rotation of water molecules. 
Due to this the proper molecular volume becomes to be dependent on the temperature and 
density. The last circumstance is proved to be essential not only for water but also for all 
liquids with non-spherical molecules and even for the noble gases, although in smaller extent. 

The comparative analysis of thermodynamic properties of the normal and heavy water 
is performed. It notices that the isotopic effect manifests itself not only in the shifts of the 
critical and melting points. In the considerably more extent it is manifested in the behavior 
of their diameters of entropy and their specific heats. 

It is taken into account that in the fluctuation region the volume, occupied by two water 
molecules, is greater than that volume, which is created by a rotating dimer. As a result 
the strong dimerization in water is stimulated. Using the methods of chemical equilibrium, 
we study the degree of dimerization. Applying the van der Waals equation of state, we 
show that the location of the critical point in water is mainly determined by the ensemble 
of dimers. It is established that different values of the specific heats for normal and heavy 
water in their fluctuation regions, 0.85 < t < 0.98, are connected with the inner rotation of 
monomers, entering the dimers {D20)2 in the heavy water. 
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II. THE DIMERIZATION OF WATER IN THE NEAR CRITICAL REGION 



In this section in accordance with [? ] we use the vdW based model with the account of 
dimerization. We calculate the average number of if-bonds per molecule along the coexis- 
tence curve. Note that for completely dimerized water hh = 1. In accordance with said in 
Introduction we assume that water near critical point is fully dimerized. Thus water dimer 
(7/20)2 plays the role of a "particle". According to the quantum chemistry calculations [? 
] the diameter of a dimer (Jd ~ 3.5A and the dipole moment dd ~ 2.6 D. 

The idea that the dimers are predominant in near critical region has been put forward 
before (see [? ] and references therein). Though the inference was made from the analysis 
of structural data (Raman scattering) which did not led to the firm conclusion. Lately, the 
results of MD simulations [? ] favored the picture that in near critical region water behaves 
more like associative fluid as opposed to the assumptions about peculiar (fractal) structure 
of hydrogen bonded network . 

Therefore we can expect that in near critical region water can be reasonably described 
by the vdW EOS for the system of dimers. The basis for the use the van der Waals like 
EOS in near critical region is due to low density and thus intensive rotation. Indeed, the 
average spacing between dimers at the critical density is about 3.5 A and is more than the 
size of dimer. However at removing of the critical point along liquid branch the character 
of the dimer rotation changes: the quasi-free rotation of a dimer discontinues at Vcav Vd- 
When Vcav{t) < Vd-, the dimers breaks and linear molecular chains are formed. A molecule 
entering in some chain can change its spatial orientation in consequence of two mechanisms: 

1) the rotation around two H-bonds connecting it with other molecules in this chain and 

2) the break of these H-bonds. passes to the rotation around the H-bond connecting dimer 
molecules. At the further removing of the critical point linear molecular chains begin to 
cooperate in spatial clusters. Here a molecule can change its orientation only because of the 
break of H-bonds connecting it with the nearest neighbors. 

The comparison of the heat of evaporation data gives the additional support for the 
picture of rotating molecular dipoles. It is known that the heat of evaporation of water is 
very big in comparison with other homologs like H2S or nonpolar simple liquids like Ar. 
This ratio for a number of the liquids is shown on Fig. |^ with H2S chosen as the reference 
fluid. Note that water and methanol have the same ratio. 
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FIG. 2: The temperature dependencies of the ratios Rnii) (a) and Rn{i) (b) according to data of 
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Note that the ratio of the heat of evaporation per molecule is almost constant for different 
fluids along the whole coexistence curve. The same constancy of the ratio of specific volume 
takes place for water and argon (see [? ]). This fact seems to be quite surprising because the 
interparticle interaction differs for different liquids. As we see, the ratio of the densities for 
water and its nearest homologue linearly increases within 10%. This increment is slightly 
larger than for corresponding ratio for water and argon (see Inroduction and [? ])• This fact 
means that the rotational motion of the weakly asymmetric molecules of liquid influences its 
EOS. This effect, as it follows from Fig. |^(b) has similar character for other liquids. Below 
it will be connected with the softness of the molecular volume entering the vdW EOS and 
close ones. 



To explain the constancy of the ratio Rnii) 



as well as Rnii) 



we are 



forced to suggest that the simplest thermodynamic quantities are determined by the averaged 
inter-particle potentials. In fact this averaged potential appear to be of the van der Waals 
character. At that time, such quantities as heat capacity, the density, the entropy diameters 
and dielectric permittivity are determined by the deviation from the averaged interactions. 
These quantities strongly depend on the short range correlations. 

To explain the numerical values of the ratio Rnii) for different liquids, we can use the 
following qualitative reasons. By order of magnitude, the heat of evaporation is the work per 
molecule which is necessary to expand its specific volume from vi to Vg-. H^'^^ ^ P {vi — Vg) , 
where all quantities are taking on the coexistence curve. Outside the fluctuation region 
vi > {~^)vg, therefore we can approximately write if*^*-* ^ Pvg or i/W ^ Z^'^tP Z^'\ where 
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TABLE I: The comparison of Rnii) calculated according to Eq. (Q) for ensemble of monomers 



(subscript (m)) and with the help of experimental data |16 |. 



Z is the dimensionless form for the compressibility factor Z = From here it follows that 
Rh is estimated as: 



rp{H2S) ^(H2S) 2{H2S) ' 

For the rare enough vapor phases of i-th and H2S system 

— ^ 1 , 



(3) 



so: 



rp{H2S) ^(H2S) 



(4) 



As is seen from Table |I| the raw estimate Eq. gives unsatisfactory prediction for liquids 
like water and methanol. Like water, the last is known of its associative properties |T^, |18 . 
Note that the compressibility factor is calculated assuming that only the monomers is present 
in near critical region. The value of Ruii) for water with account of dimerization degree is 
2.5 (see Appendix ^ and Eq. (|A4|) ) . Thus it allows to conclude that the dimerization for 
the methanol as well as water is essential in near critical region (see below). 

From our consideration it follows that the determination of the energy of H-bond for 
water from the heat of evaporation data is not correct. In accordance with our estimates, 
in the vicinity of the critical point water can be considered as the ensemble of dimers with 
small admixture of monomers. To determine the location of the critical point for water we 
apply to the van der Waals EOS: 

nT o 



p{n,T) 



an 



1 — nh 

where p is the pressure and a and h are the vdW parameters. According to they are: 

. 27r 
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TT / U{r)r'^dr. 



(6) 



In [? ? ] it had been shown that the attractive part of the dimer-dimer potential is mainly 
determined by their dipole-dipole interactions: 



UAr) = 5^%. (7) 

Therefore the coefficient a becomes to be temperature dependent: a ^ ^, and the van der 
Waals EOS transforms to Bertlo EOS: 

TiT Qj 

The coordinates of the critical point is determined by the standard conditions: 

dp p 
on ow^ 

The coordinates of the critical point for the Bertlo EOS (H) are determined by equations: 

'^'^ -36,-2vrar " V 276, " 9 a, ' ~ ^ ' ^ ~ 12-^^2 cr^ 

where A is the parameter which determines the dipole moment di^m of the dimer and 
ddim = Ae ad- We choose the value of A so that to get the best fit for the coordinates 
of the vdW critical point which are: 

Ti^) ^ 647 , Pi^) ^ 18.0 MPa , pi^^ ^ 322 kg/m^ . (10) 

The fitting gives the value A ~ 0.19 and the dipole moment of a dimer is ddim ~ 2.91^, which 
is in good agreement with the value for dimer dipole moment ddim = ^..QD [? ]. Sure total 
dimerization is an approximation but above estimates show that it is consistent. Thus the 
difference in the compressibility factor Zc is due to the deviation from vdW EOS since the 
influence of free monomers was neglected. 

Now we reflne the estimate for the location of the critical point taking into account that 
water near the critical point is the mixture of dimers and small quantity of monomers. Let 
us introduce the degree of the dimerization as following: 

ni = {l-A)no n2 = —no, (11) 

where hq = rii + 2n2 is the number density of water molecules both free and dimerized, 
n = rii + n2 is the number density of particles (both the dimers and the monomers). The 
coefficients a and b are the vdW parameters p2|: 



b = cibi + C262 , a = an + 2 CiC2ai2 + C2 022 , (12) 



where 



2tx 



CTj , ttij = TT / Uij{r)r dr, (13) 



and the concentrations of the monomers and dimers correspondingly with 

being the number density of z— th species respectively, ai are the diameters of the particles 
and Uij is the potential of the interaction between particles of i-th and j-th kinds. In [? ] it 
was shown that the interaction between monomers and dimers is determined mainly by the 
dipole-dipole forces. The corresponding potentials have the form: 

2(d?>|d|> 

«J ^ 3 ^6 

In general the virial coefficient aij is the sum of the vdW term which does not depend on the 
temperature and the temperature dependent dipole-dipole term. The chemical potential for 
monomers is: 

Hi = -T \n(l - nb) + -~2n(ci an + C2ai2 ) 

1 — no 

The chemical potential of the dimers is obtained by the change 1^2. To construct the 
binodal for the mixture of dimers and monomers we use common condition of equilibrium 
the equality of the chemical potentials and the pressures of the coexisting phases: 

l^'liq l-^vap 1 Pliq Pvap • (-^^) 



together with the chemical equilibrium equation Eq. ( [Al] ) to determine the equilibrium 
number densities of the species in coexistence phases. 

Let us discuss the compressibility factor as a function of the degree of dimerization. For 
water considered as an ensemble of monomers and described by the vdW EOS is expected 
to be 



If water is completely dimerized: 



Zi-'^w) = S/16. (16) 



This value is close to the experimental one = 0.23 near the critical point of water and 

it can be considered as one of evidences in the favor of high degree of dimerization near the 
CP. The values of Z^. corresponding to ric = Pc/ma, where is the critical mass density and 
rria is the atomic mass, for different fluids considered as the monomeric ones are in Table ||. 



As seen from Table || the compressibility factor for methanol practically coincides with 
that estimate for water which is obtained for the ensemble of dimers. This fact also points 
to the high degree of dimerization in methanol near its critical point [p3| , |2^ , that is nat- 
ural since the intermolecular interaction in methanol is determined by H-bonds, similar to 
that in water. For several other liquids in Table ^ do not forming H-bonds the values 
of the compressibility factor near their critical points are noticeably higher than for water 
and methanol. This difference is probably connected with 1) the essentially less degree of 
dimerization in them, even for Ar (see [? ]), and 2) the approximate character of the vdW 
EOS. 

Compressibility factor 
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FIG. 3: Compressibility factor as the function of dimensionless volume of the monomer 6i and 
dimensionless monomer-dimer virial coefficient ai2 (see Eq. ( [A 3D ) at the energy of dissociation 
Ed/T, = 3.5. 



The influence of the interparticle interactions on the behavior of the compressibility factor 
is presented in Fig. ^. As for the locus of the CP, the most strong dependence is observed 
on ai2. 

Ending this Section let us briefly discuss the interaction between molecules outside the 
fluctuation region. Here, as was noted above, the equation of state for water is determined by 
the averaged inter-particle potential. Besides the hard-core part it includes the contributions 
of the dispersive forces, the dipole-dipole and H-bond interactions between monomers: 

(f/(l,2)) ^ U hc{r 12) + U dis{r 12) + Umm{r 12) + UH{r 12). 
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TABLE II: The compressibility factor Z^. for fluids. 
The term Umm{j'i2) is similar to Eq. (|^): 

TJ (r\ — "'mon 

but (^^^^ < 0.1. Therefore we conclude that outside the fluctuation region f/mm(^i2) ^ 
Udisiru), while inside it Udisir^) ^ Uair). The term Uniru) plays the noticeable role only 
on small distances between monomers. In fact, it leads to some renormalization of the hard- 
core diameter. Due to this we expect that the behavior of the speciflc volume per molecule, 
the heat of evaporation and some other thermodynamic quantities outside the fluctuation 
region are determined by the potential: 

The vdW EOS (|) with the coefficient a do not depending on temperature can serve as quite 
satisfactory zero approximation for the EOS of water. 



III. THE ROTATION OF MOLECULES AND THE BEHAVIOR OF THE EN- 
TROPY DIAMETER. 

The character of the molecular motion differs in liquid and gaseous phases. The number 
density (or speciflc volume) is the natural order parameter for the liquid- vapor critical point. 
However, it is one-particle characteristics of a system and therefore it describes only crude 
enough behavior of a system. From this point of view the choice of the entropy as the order 
parameter seems to be more informative since the entropy is determined by correlation effects 
of all orders [^]. Besides the part, reducible to density, this order parameter includes also 



the contributions reflecting the more flne details of the molecular motion. Such discrepancy 
is reflected in the asymmetry of the binodal. The degree of the asymmetry of the coexistence 
curve is described by the diameter of density: 

= ^ {m + hg) - 1 , (17) 
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FIG. 4: The diameters of the entropy (a) and the density (b) for H2O and D2O according to 
experimental data |16|. 




FIG. 5: The diameters of the entropy for different hquids according to experimental data |16]. 



where hi = rii/ric ,i = I, g , ni, Ug, Uc are the number densities of the liquid, gas branches of 
the binodal and at the critical point correspondingly. The analogous quantity can be defined 
for the entropy. The temperature dependence of the entropy diameter: 



Si Sg 

Or 



(18) 



for normal and heavy waters is presented in Fig. ^. The difference in the values of Sd is 
caused by the isotopic effect (see also Section 0). The comparison of 5*^ for water with 
those for -^25', O2 and Ar, which has the spherical molecules, shows that non-monotone 
temperature dependence of the entropy diameter arises only for non-spherical molecules for 
which the rotational degrees of freedom should be taken into account. 

For all systems with non-spherical molecules there are two characteristic temperatures 
Ti and T„ at which Sd = 0. In particular, for normal and heavy waters T/ and take the 
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following values (see Fig. 



H2O : Ti = 31QK, Tu = mK, 
D2O : Ti = 280 K, T„ = 531 K . 

The nontrivial temperature dependence of Sdit) for water and substances of type H2S 
and O2 is connected with the change of the rotational motion of molecules when density of 
a system increases. In normal and heavy waters H-bonds play additional important role. 
If the average number of H-bonds per molecule is close to three or more, a molecule can 
only oscillate about temporary equilibrium position. Let tq be the average time for such 
oscillations, so called the residence time. In [? ? ? ] it was shown that the value of tq is 
essentially more than the characteristic time of the rotational motion ~ {I /ksT)^^'^ ^ where 
I is the inertia moment, only for T < Th, where Th ~ 315K. It is very surprising that the 
diameter of entropy changes its sign namely at this temperature (T; = 316i^'). At T > 
the residence time of the local molecular configurations is comparable with r^, therefore the 
rotational motion of water molecules will influence on the behavior of Sd analogously to that 
in liquid H2S or O2. 

The character of the temperature dependence becomes partly clearer when we pass to 
the excess entropy 

g{ex) ^ g _ g{id) ^ (19) 

where 

^(»^)(T) = cf)ln^-ln^, cf) = |, (20) 

1 Q Tic ^ 

is the entropy of ideal gas, k is the number of thermalized degrees of freedom which are 
treated classically. As it follows from Fig. |^ the most essential difference between the behav- 
ior of Sd for normal and heavy water as well as H2S is observed in the fluctuation region, 
where the both kind of water are close to be fully dimerized state. 

As has been said above the behavior of the diameter of the entropy reflects the change 
of the rotational motions of the molecules with the density . One can see that (heavy) water 
differs in this respect from other homologs. The situation looks like there is some deficit 
(see Fig. |^) of the order in liquid phase of water especially near the critical point. From 
the point of of the dimerization picture such a fact acquires natural explanation. The 
interaction between water molecules is converted into formation of bound states, the dimers. 
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FIG. 6: The diameter for residual (excess) entropy. 

which rotate freely in liquid phase in lower density near critical region. Note that the effect of 
lowering the entropy of the liquid phase due to nonspherical shape of the repulsive molecular 
cores was investigated by [^] within generalized vdW EOS model in low pressure region. 
Here we show that the entropy diameter is the fine characteristic which is very sensitive to 
such effect and can be described within the model EOS where the softness of the particle 
core is taken into account (see Section fV\j . To explain the behavior of the diameter of the 
entropy we use the basic thermodynamical representation of the entropy in dimensionless 
units (in units of the Boltzmann constant /c^): 

S = S, + cJn^ + f{n)~f{n,), (21) 

^ c 

where Sc is the entropy at the CP, c„ is the dimensionless specific heat, f{n) is some function 
describing the density dependence of the entropy. From Eq. (pi]) one can conclude that: 



s,.£fl±£!^,4./(-)±i(!l^-/K). (22) 



As it follows from Eq. (^2]), the behavior of the diameter essentially depends on the specific 
form of the equation of state, which determines the function f{n). In its turn, the EOS 
is sensitive to the effective proper molecular volume corresponding to the rotating non- 
spherical molecules in liquid and vapor phases. The detailed analysis of this question will 
be given in the following Section. 

Here we restrict ourselves only by the consideration of Sd{T), for which f{n) is taking 
in the ideal gas approximation. The branches of the binodal are taken in the mean field 
approximation: 

ni-ng = 2b\T\'/' + o{T), r=^^^, (23) 
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nd = a\T\+ o{t). (24) 

From here it follows that 

ni,g = l±b\T\^/^ + a\T\ + ... (25) 

Supposing that outside the fluctuation region the specific heat capacities of liquid and gas 
phases take the constant values ci^^ = | and ch^ = oif^ + Ac^,, from Eqs. ( P^ and ( P3| ) we 
get the following approximate expression for 5*^: 

Japp)^rj.^_ ( ^Cv . (id) \ ,Ti f f Ti 



Sr'iT)={ ^ + cr In^-lnM-A ( 1-^ ) ) , (26) 



As it follows from [jT6[ for Ar the value of Acy is 1, for H2S Ac^ ^ 2 and for water (both 
normal and heavy) Ac^ ^ 5. The parameter A depends on the specific EOS and 

A = 6^ - 2 a . 

For example, for the vdW EOS a = 0.4, 6 = 2, for CS EOS a = 1.2, b ^ 2.5. The behavior 
of S^^'^'^\t) for different values of A is presented in Fig. |^. Here it is necessary to note 
that Sf"'\T) becomes to be the non-monotone function of temperature for A; > 3, i.e. this 
peculiarity can be caused by both the rotational and vibration degrees of freedom. Of course, 
far away from the critical point, the contributions of the rotational degrees of freedom are 
determinative. However, already in a crude approximation described above we are able to 
estimate the position of the lesser root T/ of the equation SdiT) = 0. We see that despite the 
crudeness of the approximation Eq. (|25|) for the binodal at suitable values of A the reasonable 
values for the root Ti are obtained. 



In the Section it will be shown that the behavior of Sd{T) is noticeably changed if 
f{n) is taking in the Carnahan-Starling form with modified proper molecular volume. 

We want to finish this Section by the brief discussion of physical cause which leads to the 
appearance of the upper root for the equation Sd{T) = 0. The definition Eq. (0) of the 
entropy diameter includes the value of entropy Sc at the critical point. It is known that the 
value of Sc"^^^ calculated in the mean field approximation (for example in the vdW or CS 
approximations) is higher than the experimental value S'f^^\ It is a result of the neglecting 
by the mode-mode interaction in the mean field approximation. Therefore for the correct 
comparison of Sd{T) calculated on the basis of experimental data with its mean field analog 
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(a) (b) 

FIG. 7: S^f^^^ according to Eq. (|2|) as a function of A (a), (b) the comparison of Eq. ( p6[ ) (solid 
line) with the parameters Ac^ = 5 and A; = 6, a ~ 0.73 , b ~ 2.3 , A ~ 3 corresponding for water 
with the data (dots) for Sd [|^]. 

s^J^^\t) it is necessary to pass to the modified entropy diameter Sjf^^ = S^'^^ +ASc, where 
ASc = S'c™'^^ — Sc- It is clear that the equation S''/'^\t) = will take two roots: the value of 
Ti is shifted and the additional root T„ appears. This situation will be considered in more 
details in separated work. 

IV. THE "SOFT" CORE EOS MODEL 

In this Section we discuss the nature of the non-monotone temperature dependence of 
the diameter of entropy. We will show that the region of negative values of Sd is mainly 
caused by the variation of the molecular volume in dependence of pressure. The appearance 
of such dependence is natural for liquids with non-spherical molecules. Indeed, the proper 
molecular volume for them is determined as the volume of cavity formed by arbitrary rotating 
molecule. If the rotation of molecules becomes restricted, the corresponding cavity volume 
diminishes. Obviously, such an effect should take place with increasing of density. In 
water the considerable additional influence on the value of the effective molecular volume is 
generated by H-bonds. To describe the influence of the effective molecular volume on the 
behavior of the diameter of entropy, we start from some EOS, in which the dependence of 
the molecular volume on pressure is taken into account. Then we reconstruct the free energy 
and after this we find entropy in vapor and liquid phases and also the diameter of entropy. 
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The EOS for simple molecular fluids is usually taken in the form p2| ? ]: 

p = p+in,T) +p_{n,T) , 



(27) 



where p^ is the pressure contribution due to hard core repulsive interactions and p_ is the 
contribution of the attractive long range part of the potential. For simple fluids it is taken 
in the standard form: 

p_ = -an^. (28) 

From what has been said above it is clear that one can model the influence of the density on 
the rotation as the restriction of the angular conflguration space available for the molecule 
in a cage formed by its neighbors. The change in the available space for free rotation 
of the molecule can be described by the dependence of the volume parameter b on the 
density and the temperature. In accordance with said above we suppose that the density 
and temperature dependencies of the effective molecular volume b can be modeled by the 
relation: 

^° p(*^)=nf. (29) 



This parameter is included in the corresponding pressure term which describes the short 
range interaction (the hard core). The modifled vdW EOS approximation for the free energy 
due to the interaction: 

p = = = + riT 30) 

(1 + 7nT)(3(l + 7nT) -n) 8 

For small 7 the coordinates of the critical point are: 

Tcil) 19 , , rich) , , 

The analogous relations for EOS corresponding to the following model of the free energy: 

inQ\ ^bn ( A — bn) , , 

F^^^^ = F,d + T ^ /-an, 31 

where Fid is the ideal gas term are: 

^ = l + t,7 + o(7), !^ = l + ni7 + o(7), ti ^ 4.5 , m ^ 8.1 

The modiflcation of the binodal for the Carnahan-Starling EOS at 7 7^ is shown on 
Fig. ^(b)| . As we see that the width of the binodal decreases with 7. This naturally reflects 
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FIG. 8: The dependence of the coordinates of the CP (a) and the binodal (b) on 7 for the model 
Eq. (|lD with the soft core Eq. (p|). 



the diminishing of the phase asymmetry due to bigger compressibility of the effective volume, 
which in its turn means that the asymmetry between particle and the cavity (hole) decreases. 
From such point of view the most symmetrical is the binodal for liquid helium [? ]. The 
calculation of the entropy diameter 5*^ supports this conclusion (see Fig. 

Entropy of the system, corresponding to the Carnahan-Starling EOS, is determined in 
the standard way: 

S = -^-4^ (32) 



dT 



The corresponding diameter of the entropy is presented in Fig. |ll(b)| . Note that for water 
we neglected the temperature dependence of the coefficient a in Eq. (pSD in accordance with 
the discussion in Section |Tl|. Thus the EOS with the soft core forms the adequate basis for 
the successful description of the binodal and fine nontrivial details in the behavior of the 
entropy diameter. 



V. MANIFESTATION OF ISOTOPIC EFFECT 



In this Section we want to complete our results presented above by the brief discussion 
of some distinctions in the behavior of the normal and heavy water. We focus our attention 
on the following characteristic manifestation of the isotopic effect: 1) small differences in 
the locations of the triple and critical points; 2) essentially more noticeable differences in 
the locations of the roots ti and tu of equations Sd{t) = (t^ = Tk/Tjt\ k = l,u) and 3) the 
nontrivial difference in the behavior of the specific heats. In these cases we observe different 
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FIG. 9: The dependence of on 7 for Eq. (|3l[). 
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FIG. 10: The calculated binodal (solid) for the model Eq. ( pl| ) with 7 k, 0.03 and the binodal for 
R2O and H^S (dots) 
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FIG. 11: The entropy diameter 5^ for the model Eq. (31) with 7 ~ 0.03 (solid line) and the data 
for and i?2'S' (dots). The shifted entopy diameter curve (see Section |T|) is shown by the 

dashed line. 
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manifestations of the isotopic effect. 

The critical temperature of normal water, Tc{H20) = 647K, is rather higher than one, 
Tc{D20) = 644K, for heavy water. It means, that the position of the critical point in heavy 
water is determined by the ensemble of dimers doped by small quantity of monomers, i.e. 
similarly to that for normal water. The small lowering of the critical temperature for heavy 
water is probably connected with 1) slightly different values of parameter A^/cr^ for H2O and 
1^20 (see Eq. (P) and 2) some distinction of the dispersive interactions between molecules. 
Note, that in Eq. (P) the contribution of the dispersive forces is at all ignored. 

The difference in the locations of the triple points is also insignificant: Ttr{H20) = 273K 
and Ttr{D20) = 277K. From here it follows that normal and heavy water belong to the same 
class of corresponding states. For them Ttr{H20) /T^{H20) = 0.42, Ttr{D20)/T^{D20) = 
0.43, that essentially differs from the ratios for the water homologue H2S and benzene: 
Ttr{H2S) /Tc{H2S) = 0.50, TtriC^HQ) /Tc{CqHq) = 0.50. Such significant discrepancy is con- 
nected with the considerable influence of the dimerization on the location of the critical 
point. From Fig. |l| it follows that due to dimerization the specific volume per molecule 
increases approximately on 15% . Assuming that the temperature shift has the analo- 
gous value: ~ we obtain the following estimate for the critical temperature of 
non-dimerized water (Tc""^^ ^ — AT^): Tc^'^\H20) = 550K. It leads to the estimate: 
Ttr{H20)/T^'"^\H20) ^ 0.50. This fact can be considered as an indirect evidence for the 
dimerization near the critical point. 

The locations of the roots ti and tu of equations Sd{t) = for normal and heavy water 
satisfy the following relations: 

ti{H20) - ti{D20) ^ 0.053, 

K{D20) - K{H20) ^ 0.054. 

Their surprising closeness is not occasional. In the Section 4 it was stressed that the non- 
monotone temperature dependencies of the entropy diameter S'^\t), k = H20,D20, is 
connected with the rotational motion of molecules. Therefore the comparison of the dif- 
ferences |Ati| ^ |At„| ^ 0.053 with the characteristic parameter /i = (-^^o) -a;^(£'2 o) 
the rotational motion is quite relevant. In principle, molecules can change their relative 
orientations by many ways; 1) each molecule rotates independently of others; 2) two nearest 
molecules change their orientations in concord, while all other neighbors retain in initial po- 
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FIG. 12: The relative value of the difference of the specific heats for normal and heavy water 



according to |16]. 



sitions and so on. Note, that the dimerization in vapor phase of water leads to the natural 
lowering of entropy and negative values of Sd{t). In liquid phase the existence of strong ori- 
entation correlations makes improbable the first way. At the same time, the corresponding 
pair in the second case can be identified with a dimer, for which we find an estimate: 



/X ~ 1 — ./ ^ 0.05. 

V ITT'D^O 

It is necessary to emphasize that we do not speak about the rotation of pairs (dimers) 
in liquid phase. Hence, the very close values of two dimensionless parameters is a strong 
argument in the favor of the second scenario of reorientations in vapor and liquid phases of 
normal and heavy water. 

It seems to be natural to suppose that similar character of the rotational motion in normal 
and heavy water will also lead to the following estimate for the ratio of their specific heats: 



This ratio is really in correspondence with experimental data, presented in Fig. |T2| However, 
here we would like to focus our attention on the practically constant value of the difference 
Ci^^^\t) — Cv^^^\t) in the temperature interval 0.8 < t < 0.95. It corresponds to the 
fluctuation region, in which the normal and heavy waters are strongly dimerized. As we 
see, in this region Ci^^^\t) — Ci^'^'^\t) ^ i. This effect finds the natural explanation if we 
assume that the inner rotation in dimers {D20)2 is possible (see Fig. |^). The activation 
of this type of thermal motion in heavy water is justified by greater length of D-bonds in 
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FIG. 13: The relative rotation of the monomers entering a dimer. The parameters of water dimer 



are taken from wwwl. lsbu.ac.uk/water/ 
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FIG. 14: Temperature dependence of the ratio Rv{t) 
curves of the normal and heavy water ||l6||. 



0.8 



0.9 



{t) / v^^'^^\t) on the coexistence 



comparison with H-bonds |^ and low energetic barriers ^3|. This circumstance becomes 
evident from the comparison of the specific volumes for normal and heavy waters (Fig. 
The values of Rvit) near the triple points, Ur ~ 0.42, are mainly caused by different lengths 
of H- and D-bonds, since the numbers of them per molecule in normal and heavy water are 



practically the same [? ] . From Fig. |14| it follows that the length of a D-bond is greater than 
one for an H-bond approximately 1.5%. 



VI. CONCLUSION 

The main attention in this paper is focused on the physical nature of the temperature 
dependence of main thermodynamic properties of water: density, heat of evaporation and 
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entropy as well as specific heat everywhere on the coexistence curve. All temperature interval 
of the two-phase vapor-liquid states is considered. 

Analyzing the behavior of the density and the heat of evaporation of water we have estab- 
lished that these quantities have argon-like behavior everywhere excepting the fluctuation 
region near the critical point. It means that such a character of the temperature depen- 
dence is formed by averaged argon-like inter-particle potential, which arises in consequence 
of the rotation of water molecules. H-bonds bring in the specificity in the rotational motion 
of molecules, however their influence is not crucial. This circumstance is justified by that 
fact that the shear viscosities of water and argon have the same order of magnitude. The 
situation is considerably changed only for supercooled states of water and near its critical 
point. In the present paper wc touch some details only for the near critical behavior. 

Here the volume occupied by two water molecules becomes greater than that volume, 
which corresponds to a rotating dimer. As a result, the favorable conditions for the dimer- 
ization are created. Applying the method of chemical equilibrium it is shown that the degree 
of dimerization exceeds 0.9. Due to dimerization the characteristic ratio Rv{t) for specific 
volumes of water and argon increments approximately 15% in the vicinity of their critical 
points. It shown that the ensemble of dimcrs with small admixture of monomers allows to 
reproduce successfully the location of the critical point of water. It is also responsible for 
the correct relations between critical amplitudes. 

At elongation from the fluctuation region dimers destroy and short-living linear chains of 
water molecules are formed. Here the thermodynamic properties of water are determined 
by rotating monomers (of course, the rotation of monomers is not free). Due to formation 
of short living H-bonds the effective proper volume, corresponding a monomer, becomes to 
be dependent on the temperature and density, or pressure. Including the dependence of 
the proper molecular volume on the pressure to the van der Waals or Carnahan-Starling 
equations of state leads to the very nontrivial consequences. Among them we mark out 1) 
the non-monotone temperature dependence of the diameter of entropy and 2) the differences 
in the values of the specific heats for vapor and liquid phases, calculated with the help of 
the modified van der Waals or Carnahan-Starling equations of states. Here it is appropriate 
to note that the specific heats for coexisting phases arc identical if the proper molecular 
volume is constant. Besides, the change of the sign of the entropy diameter at the lower 
limit for the fluctuation region tu ~ 0.8 is not also occasional. It is a strong argument in the 
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favor that this root of the entropy diameter is connected with the discrepancy of the mean 
field and experimental values of entropy at the critical point. 

The diameter of entropy changes also its sign at the lower temperature t; ~ 0.5. The 
physical prerequisite for this root of the entropy diameter is the following: near ti the ensem- 
ble of linear molecular chains begins to form the spatially arranged the H-bond network. The 
rotation of molecules becomes more difficult and the character of the density dependence of 
the proper molecular volume essentially changes. 

The very nontrivial manifestation of the isotopic effects is observed at the comparison of 
the temperature dependencies for the diameter of entropy and the specific heat of the normal 
and heavy water. Here we pay attention on the characteristic values of differences in the 
positions of the upper and lower dimensionless temperature for the normal and heavy water. 
These differences are close to each other and they practically coincide with dimensionless 
parameter describing the discrepancy in the rotations of dimers arising in vapor and liquid 
phases of normal and heavy water. Another remarkable fact is the difference of the specific 
heats for the normal and heavy waters in the fluctuation temperature interval 0.85 < t < 
0.98. In accordance with our reasons it arises owing to switching on the internal rotation 
for dimers in the heavy water. 

The small differences in the values of the crystallization and critical temperatures for 
them is connected with some distinction of their averaged intermolecular potentials. In 
consequence of the strong dimerization near the critical point the normal and heavy water 
form the separated class from the point of view of the principle of corresponding states. 
Besides water, many alcohols belong probably to this class of corresponding states. 

APPENDIX A: THE DEGREE OF DIMERIZATION OF WATER NEAR THE 
CRITICAL REGION (MEAN FIELD ANALYSIS) 

The equation for the chemical equilibrium between dimers and monomers is: 



where fi^ are the excess chemical potentials (with respect to the ideal gas term) of the 
monomers and dimers correspondingly which are determined from the excessive part of the 
free energy corresponding to the EOS (^. The quantity K[T) is the constant of chemical 




(Al) 
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equilibrium we can approximate it by simple expression [p^: 



K{T) =exp(^^y (A2) 

The dissociation energy Ed corresponds to the energy of H-bond which in kBTc units is 
of order Eh = 3 -t- iksTc |]3^, Such a large value of the H-bond energy allows us to 
conclude that near critical water is dimerized to a great extent [? ? ]. In particular it agrees 



with the results of the numerical simulations 37 



In order to calculate the dimerization degree at the critical point we use the model Eq. (H) 
with the dissociation considered as the perturbation to the completely dimerized state. In 
this way we obtain the critical temperature in the linear approximation on e~^'^/'^ < 1: 



1 2 



where 



The corresponding value for Zc is: 

Z, = zW(l-A,/2) . (A4) 

The dependence of Z^. on the relevant parameters is shown on Fig. |^. E.g. at ai2 ~ 0.2, 
bi = 0.2 and Ed/Tc = 3.5 it takes the value Zc ~ 0.2. In the vicinity of the critical point 
assuming that e <^ 1 , from (|AT|) we get the degree of the dimerization: 



A, = 1 - - e-^''/2-| + O (e^) , (A5) 

3 

where = j is the parameter which determines the difference in charge distribution 

in the monomer and the dimer. The dependence of the dimerization factor on the parameters 



of the monomer-dimerization interaction is shown on Fig. |15. 

Thus in the critical region the degree of dimerization A at ai2 ~ 0.4, which corresponds 
to di/d2 ~ 0.5,o"i/cr2 ~ 0.7, is approximately 0.9. This conforms with the thermodynamic 
estimates given in Section | and the results of works [? ] . 
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Dimerization degree 




FIG. 15: Dimerization degree as the function of the parameters of monomer-dimerization interac- 
tion 
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